Did supermassive black holes form by direct collapse? 
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Abstract. Rapid infall of gas in the nuclei of galaxies could lead to the formation of black holes by direct collapse, without 
first forming stars. Black holes formed in this way would have initial masses of a few M s , but would be embedded in massive 
envelopes that would allow them to grow at a highly super-Eddington rate. Thus, seed black holes as large as 10 3 — 10 4 Mq 
could form very rapidly. I will sketch the basic physics of the direct collapse process and the properties of the accreting 
envelopes. 
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INTRODUCTION 

Despite years of study, we still do not know how the 
seeds of supermassive black holes formed. Few if any 
of the pathways in Martin Rees's famous flow chart 
(Begelman & Rees 1978) can be ruled out, but none of 
the routes is particularly well understood, either. What 
we do know is that some very massive (> 1O 9 M ) black 
holes had to exist by z ~ 6 in order to explain early 
quasars (Fan 2006). If the seeds of these black holes were 
the remnants of massive stars, then they must have grown 
by Eddington-limited accretion for most of the time since 
their formation, or else much of their growth was due to 
mergers. A second possibility is that the seeds formed 
by such a rapid accumulation of matter that it may be 
considered to be a direct collapse. I will focus on the 
latter possibility in this paper. 

DIRECT COLLAPSE 

The Pop III star formation processes we have heard 
about at this conference result from the infall of gas 
at rates ~ 10~ 4 — 10~ 2 M o yr _1 . What would happen 
if the infall rate were much higher? The entropy of 
matter laid down by gravitational infall onto a growing 
central mass increases with time. At low inflow rates, 
however, nuclear ignition halts the contraction of the 
core and raises the entropy in the interior, leading to a 
high-entropy object — a star. If the inflow rate is high 
enough, however, the core will be so tightly bound by 
the time nuclear reactions start that the energy release 
will be insufficient to halt core contraction. In this case, 
we are left with an object with a low-entropy core and 
a high-entropy envelope. This is the situation that can 
lead to the direct formation of a black hole, without a 
stellar precursor. This situation should apply when the 
infall rate exceeds a few tenths of a solar mass per year 



(Begelman et al. 2006), although more work needs to be 
done to refine this estimate. 



Inflow rate 

The conditions under which such high inflow rates 
might occur are very uncertain. Such rapid infall would 
almost certainly be driven by gravitational torques, 
which could be local (Gammie 2001) or global (as in the 
"bars within bars" mechanism: Shlosman et al. 1989). 
The "natural" gravitational inflow rate is given by M ~ 
V 3 /G ~ 0.1(v/10 km s ) 3 M Q yr _1 , where v represents 
the internal velocity dispersion (~ turbulent or sound 
speed) for a locally unstable thin disk and the orbital 
speed for a globally unstable (fully self-gravitating) sys- 
tem. This means that very large inflow rates are possi- 
ble in dark matter haloes with velocity dispersions ex- 
ceeding about 10 km s~ ! , which have masses exceed- 
ing ~ 10 8 — 10 9 M a and become common at redshifts 
~ 10 — 15. Global gravitational instabilities could occur 
in a significant fraction of such haloes if the gas is un- 
able to cool much below the virial temperature (~ 10 4 
K), which requires both that they have been spared sig- 
nificant metal enrichment and that H2 formation is sup- 
pressed (Bromm & Loeb 2003; Begelman et al. 2006). 
However, recent calculations by Wise & Abel (2007, and 
these proceedings) suggest that molecular hydrogen for- 
mation is inevitable. This may have the additional effect 
that many more haloes formed Pop III stars at higher red- 
shifts, making it harder to avoid metal enrichment and 
possibly depleting the supply of gas available for infall 
as haloes merge. 

Rapid infall may be easier in more massive haloes. 
There will be a larger disparity between the thermal tem- 
perature of the gas and the virial temperature, but pre- 
sumably the gas will form a multiphase structure with 



turbulent velocities (perhaps driven by stellar energy 
sources) dominating the internal energy — much like the 
interstellar medium. Simulations are beginning to eluci- 
date these structures, but we have a long way to go before 
we understand the tradeoffs between large scale inflow 
and in situ fragmentation and star formation. Inflow may 
be stimulated by the large-scale gravitational torques as- 
sociated with mergers, a key element in the simulations 
of black-hole growth and fuelling (e.g., Di Matteo et 
al. 2005). Finally, we note that such inflow rates must 
be possible, because they are required to power quasars. 
The open question is whether they can also occur when 
the black hole is absent, or very small. For the remainder 
of this article we will assume that the answer is affirma- 
tive. 



Black hole formation 

Once the mass of accumulated gas exceeds a few solar 
masses, radiation pressure dominates the envelope. The 
core, with gas pressure comparable to radiation pressure, 
maintains a roughly constant mass <~ 10M©. The bound- 
ary between dynamic infall and quasistatic contraction 
is close to the radius where infalling gas liberates en- 
ergy at the Eddington limit. For a constant infall rate of 
M = 0.1m_iM e yr _1 , this radius turns out to be constant, 
R em ~ lm_i AU. The core shrinks under the increas- 
ing pressure of the envelope, R cole ~ R env /m em , where 
m em ~ 0.1m_if yr is the envelope mass in solar units. 

One can sketch the likely interior structure of the en- 
velope using scaling arguments. Since the interior is 
undergoing Kelvin-Helmholtz contraction, the density 
profile adjusts so that the diffusion timescale fdiff ~ 
(pKr 2 1 'c)(T I V|Vr|) is of order the elapsed time at all 
radii. This means that the specific entropy of any mass 
shell, s(M), declines slowly (logarithmically) with time, 
and is never far below the value it had when the mass 
shell was added. In the radiation pressure-dominated en- 
velope, s(M) °< Pr/Pg x M(r) x l 2 , where p r °< T 4 and 
Pg°< pT are the radiation and gas pressure, respectively. 
There are then two possible scaling laws for the en- 
velope structure. If T/r\VT\ ~ 0(1), then p(r) °< r~ 2 , 
T(r) °< r -1 / 2 (for uniform opacity k), and the enve- 
lope joins smoothly onto the core. This was the result 
presented in Begelman et al. (2006). The other possi- 
ble scaling law has T/r\VT\ 1, i.e., nearly isother- 
mal structure. The density profile in the envelope is then 
p(r) °c r -1 / 2 , and there is a large jump in density and 
temperature going from the envelope to the core. The lat- 
ter structure seems to be the one that emerges naturally 
from "nearly self-similar" models for the core+envelope, 
which neglect rotation. Given the manner in which the 
mass accumulates, though, rotation is bound to be im- 



portant. Moreover, even a small amount of it can affect 
the structure dramatically, increasing the binding energy 
of the gas (which is otherwise very weakly bound since 
the mean equation of state has y « 4/3) and protecting 
it from pulsational instabilities. Including the effects of 
rotation will require the modeling of angular momentum 
transport and dissipation, but these complications are not 
prohibitive. It will be interesting to see which of the scal- 
ings applies when rotation is included. 

Fortunately, the behavior of the core is insensitive to 
the structure of the envelope. Since the core is largely 
dominated by gas pressure, its temperature must track 
the virial temperature, r core °= m env . Rapid nuclear burn- 
ing starts in the core when the envelope mass is ~ 100M© 
and the density is not that different from the cores of 
main sequence stars, but the system passes through this 
phase so quickly (f evo i ~ 10 3 mZ l yr) that the energy re- 
lease has little effect. By the time the envelope mass 
reaches several thousand solar masses, the core tem- 
perature has climbed to a few xlO 8 K. Because of the 
extremely deep potential well created by the envelope, 
the energy released by nuclear reactions at this stage is 
unable to unbind the gas (especially if rotation has in- 
creased its binding energy) and the core proceeds to the 
temperatures (~ 5 x 10 8 K) at which runaway neutrino 
losses occur. At this point the core loses pressure support 
and collapses to form a ~ 10 — 2OM black hole. 



POST-BH EVOLUTION 
Super-Eddington accretion 

At the time of its formation, the black hole is embed- 
ded in an envelope of more than a hundred times its mass. 
If it were limited to accreting at the Eddington limit for 
the mass of the black hole, the pressure of the escap- 
ing radiation would have essentially no effect on the en- 
velope. In fact, the pressure of the envelope should be 
able to drive the accretion rate up to the point where 
the liberated luminosity approaches the Eddington limit 
for the mass of the envelope. Relative to the mass of the 
black hole, the accretion is super-Eddington by a factor 
M em /M BH . 

Accretion inside the massive envelope (which itself 
continues to grow at a rate M) can lead to very rapid 
growth of the black hole. The Eddington ratio for black 
hole growth is 
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where 0.1e_i is the accretion efficiency. For an initial 
black hole mass of 10M o , the e-folding time is a thou- 
sandth of the Salpeter time, or only ~ 10 4 yr, and super- 




FIGURE 1. Schematic illustration of the quasistar, from 
Begelman et al. (2007). A seed black hole of mass Mbh ac- 
cretes gas from a massive, radiation pressure-supported enve- 
lope at a rate set by conditions outside the Bondi radius. The lu- 
minosity liberated by the accretion process is transported con- 
vectively in the inner regions of the envelope, with a transition 
to a radiative zone where convection becomes inefficient. We 
illustrate rotational flattening and ongoing disk accretion at a 
fraction of a Solar mass per year. 



Eddington growth appears possible up to black hole 
masses > 10 6 M e . However, this estimate does not take 
fully into account the back reaction of black hole growth 
on the envelope. We will see below that the era of rapid 
growth is limited to much smaller (but still interesting!) 
black hole masses. 



Quasistars 

The energy liberated by accretion has to escape. Since 
the accretion flow is rotating, some of it could exhaust 
through a low-density funnel. But much of it presum- 
ably percolates through the accretion flow and envelope. 
Thus, the accreting black hole provides an energy source 
for the envelope, which is therefore a kind of star-like 
object which we have dubbed a "quasistar" (Begelman 
et al. 2006). Since there is a limit to the outward energy 
flux that the accreting gas can carry, the accretion rate is 
regulated at a fraction ~ e~ l {c s /c) 2 of the Bondi rate, 
where c s <~ {p r /p) 1 ^ 2 is the sound speed at the Bondi 
radius (Gruzinov 1998; Blandford & Begelman 1999; 
Narayan et al. 2000; Quataert & Gruzinov 2000). The ac- 
cretion rate, and associated energy flux, can be expressed 
in terms of the temperature and density deep within the 
quasistar, but outside the black hole's sphere of influence. 
The energy flux must also equal the Eddington limit for 
the entire quasistar, Le{M*) (where we now use to 
denote the mass of the quasistar) — the interior structure 
adjusts so that this is satisfied. Since the Eddington limit 
is lower inside the quasistar (where the enclosed mass is 
lower), the quasistar's interior is strongly convective and 
can be modeled by an n = 3 polytrope. (Strictly speaking, 
one should use a "loaded polytrope" model [Huntley & 
Saslaw 1975], taking into account the mass of the black 
hole, but the loaded and unloaded models converge in the 
region of interest when M % ^> Mbh-) The interior struc- 
ture of a quasistar is shown schematically in Fig. 1. 



According to the polytropic relations, the density and 
temperature outside the Bondi radius are uniquely related 
to the mass and radius (/?*) of the quasistar. Expressing 
the accretion rate in terms of these quantities and setting 
the luminosity equal to Le{M % ), we can express R t in 
terms of M* and Mbh- 

fi*~4x lO l4 a 2/5 m 4 J^m^ l/5 cm, (2) 

where the masses are expressed in solar units and a < 1 
parametrizes the efficiency of angular momentum and 
energy transport inside the Bondi radius (Begelman et 
al. 2007). The photospheric temperature is 

r ph ~ 10 3 a- l/5 m B H /5 ml /2 ° K. (3) 

Since the black hole growth rate is proportional to M*, 
while the quasistar mass increases linearly with time, 
Mbh 00 M 2 at late times. This implies that the quasis- 
tar's radius grows with time, and its photosphere be- 
comes cooler. (The result holds even if the quasistar mass 
is fixed.) The appearance of a quasistar differs dramati- 
cally from that of the pre-black hole envelope. Shortly 
after the black hole forms the envelope expands to ~ 100 
AU (from ~ 1 AU), and the interior temperature drops 
to ~ 10 6 K, quenching all nuclear reactions. A quasis- 
tar resembles a red supergiant, except that it is radia- 
tion pressure supported and its energy source is accre- 
tion. In some respects it is reminiscent of a very mas- 
sive Thorne-Zytkow (1977) object, but with crucial dif- 
ferences. Besides being dominated by radiation pressure 
and powered exclusively by accretion, quasistars have a 
distributed rather than shell-like energy source, which 
cannot be regulated by the slow settling that character- 
izes accretion onto a neutron star. As a result, quasistars 
come to grief if their photospheric temperatures get too 
low. 

Begelman et al. (2007) computed quasistar models us- 
ing Pop III opacities from Mayer & Duschl (2005). As 
in red giants with standard abundances, there is a min- 
imum photospheric temperature associated with a sharp 
drop in opacity. Because of the absence of metals, the 
"opacity crisis" occurs at a somewhat higher temperature 
— around 4000 K — than for metal-enriched compo- 
sitions. Like the Hayashi track (Hayashi 1961; Hayashi 
& Hoshi 1961) for red giants and convective protostars, 
the minimum temperature of quasistars arises from the 
impossibility of matching the convective interior to the 
radiative zone and photosphere. The details are some- 
what different, however, because radiation pressure dom- 
inates in quasistars, whereas the convective zones of or- 
dinary red giants are gas pressure-dominated and resem- 
ble n = 3/2 polytropes. As a quasistar crosses into the 
forbidden zone, the flux escaping from the convective in- 
terior exceeds the Eddington limit and prevents the qua- 
sistar from maintaining hydrostatic equilibrium. Fig. 2 
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FIGURE 2. The shaded areas show forbidden photospheric temperatures of Pop III quasistars, as a function of envelope and 
black hole mass, from Begelman et al. (2007). The lighter shaded region is computed assuming that a = 0.1, the darker shaded 
region is for a = 0.05, while the dashed line is for an analytic "toy" opacity model. Superimposed on the figure are evolutionary 
tracks for an accretion rate onto the envelope of M„= 1 M© yr _1 and M*= 0.1 M yr , The evolution pushes the quasistar into 
the forbidden region of parameter space, where it evaporates. 



shows the forbidden zone in the M* — Mbh plane, along 
with representative evolutionary tracks. 

Unlike ordinary red giants reaching the Hayashi track, 
quasistar photospheres cannot hover stably at close to the 
minimum temperature. Fixing the photospheric temper- 
ature, we find that the Eddington ratio scales as 



(4) 



Thus, growth of the black hole exacerbates the dynam- 
ical imbalance, as does partial evaporation of the enve- 
lop43 — the opacity crisis is an unstable situation. We 
predict that quasistars entering the forbidden zone must 
evaporate rather quickly. 

There is a complication due to a bump in the opacity 
from bound-free transitions, which creates a narrow re- 



1 Since the black hole's growth rate is proportional to M t , rapid accre- 
tion onto the envelope also would not be able to stabilize the quasistar 
for very long. 



gion in the radiative zone where the luminosity is super- 
Eddington. A similar problem arises in models of Lumi- 
nous Blue Variables (e.g., Owocki et al. 2004), and re- 
mains unresolved. However, it does not seem likely that 
this feature alone will lead to catastrophic mass loss in 
the case of quasistars. 



DISCUSSION 

The rapid infall of gas in galactic nuclei or pregalactic 
haloes provides a means for forming seed black holes 
and rapidly growing them into the intermediate mass 
regime. For Pop III quasistars with inflow rates ~ 0.1 — 
IMq yr~', and simple assumptions about parameters 
like a, the black holes could reach masses ~ 10 3 — 
1O 4 M before the quasistar evaporates. The quasistar 
masses could be as large as 10 5 — 1O 6 M0. Metal-rich 
quasistars could reach somewhat lower temperatures, but 
it is not clear whether this implies larger black hole 



masses, because the run of opacity in the radiative zone 
would be more complex and dust formation might lead 
to enhanced mass loss. 

Are quasistars detectable? If so, they would be seen 
at their most massive, shortly before they evaporate. 
Spectrally, they would resemble 2000-4000 K blackbod- 
ies, depending on metallicity, and in the Pop III limit 
their spectra would be featureless. Since they radiate 
at the Eddington limit corresponding to the opacity at 
the convective-radiative transition — which is close to 
that of electron scattering — their luminosities could 
reach 10 43 — 10 44 erg s _1 . However, their short lifetimes 
(~ 10 5 — 10 6 yr) would make them fairly rare. 

Quasistars can exist only when the envelope mass 
greatly exceeds the mass of the embedded black hole. 
Therefore, they are unlikely to form in the nuclei of 
galaxies that already possess a supermassive black hole 
— in this case, a period of rapid infall (e.g., following a 
merger) would presumably trigger a quasar outburst in- 
stead. However, low-redshift quasistars might conceiv- 
ably form in galactic nuclei in which the black hole had 
been ejected due to three-body interactions, or in which 
a black hole had never formed. The earliest plausible 
sites of quasistar formation would have been pregalactic 
haloes with virial temperatures exceeding 10 4 K. These 
would have been most common at redshifts ^6 15. 
The spectra of Pop III quasistars at these redshifts would 
peak at about 10/xm, but they might be marginally de- 
tectable by James Webb Space Telescope on the Wien tail 
at 3 — 5/im. If direct collapse in 10 4 K haloes is a prin- 
cipal route for forming supermassive black hole seeds, 
there could be as many as 1-10 per JWST field, but iden- 
tifying them would be an extreme challenge. 
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